Massive Free-Streaming Neutrinos and Rise of N v at Recombination 
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We present the Einstein- Vlasov solution for the momentum distribution of the relic free-streaming 
neutrinos. We show that it is possible to explain the rise in the effective number of neutrinos N v 
from those present at the end of the big bang nucleosynthesis (BBN) N„(Tbbn) = 3.7lig.4s towards 
N v (T r ) = 4.34l °gg noted at time of electron-ion recombination (r). This increase is due to the 
interplay of the neutrino mass and the non-equilibrium form of the neutrino distribution after the 
freeze-out. One implication of our scenario is that the mass of the heaviest neutrino should be 
near the recombination temperature, T r = 0.253 eV. If instead of the Einstein- Vlasov solution, a 
thermal equilibrium distribution is inadvertently invoked, one would expect a decrease in N v at 
recombination. 

PACS numbers: 51.10.+y,95.30.Cq,14.60.Pq,26.35.+c 



I. INTRODUCTION 

The relic (i.e. background cosmic) neutrinos 
have not been directly measured [l], [|| . Their pres- 
ence and properties are inferred from reaction dy- 
namics throughout the history of the universe [!, Q • 
Chemical freeze-out occurs when particle number 
changing reactions stop due to the expansion of the 
Universe and a decrease in particle densities near 
T c h — 3.3 MeV. All 2-into-2 particle collisions involv- 
ing neutrinos cease near Tf ~ 2.3 MeV and neutri- 
nos begin to free-stream, akin to the free-streaming 
of photons after atoms form and ion-electron plasma 
entirely disappears at the much lower recombination 
temperature T r = 0.253±0.001 eV Hi]. In a future 
work we will present a detailed study of the differ- 
ent chemical and kinetic freeze-out neutrino condi- 
tions 0]. 

After neutrino abundance i.e. chemical equilib- 
rium condition is reached (chemical freeze-out) there 
are some known processes that feed into the neutrino 
gas and result in a post-equilibrium addition to the 



neutrino abundance, an effect that has been stud- 
ied extensively in the literature [§| . To quantify this 
effect one introduces effective number abundance of 
neutrinos of N u defined as the relic neutrino to the 
photon energy density ratio, allowing for photon to 
massless fermion statistical ratio 7/8, to be com- 
pared to the standard model which contains exactly 
N v = 3 left handed neutrino flavors. 

The modification due to two body reactions 
involving neutrinos and from neutrino chem- 
ical (abundance) freeze-out until the time when 
positrons have annihilated was calculated in Q and 
we will often refer to this value N^(Tbbn) = 
N thBBN _ 3 046. Neutrino feeding occurs in the 
BBN phase finishing at Tbbn — m e /10 but the 
net result is smaller [l0| . All three flavors are 
strongly mixing as a function of time. In the 
early universe environment as long as at least one 
of the neutrinos remains in equilibrium with the 
background, all three flavors participate in the ki- 
netic processes equilibrating with the background, 
consisting mainly of photons and electron-positron 
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plasma. 

There is currently a significant degree of in- 
terest in the precise value of N v . Studies of the 
cosmic microwave background (CMB) and element 
production during big bang nucleosynthesis (BBN) 
put constraints on the effective number of neutri- 
nos. In the WMAP 7-year report [ll| the effective 
number of neutrinos at recombination T r — 0.253 
eV is N v (T r ) ee Nl = 4.34± ^, while the current 
value deduced from BBN is NJT BBN ) = N BBN = 
3.71^^45 [3,[l3l- The errors inherent in the analyzed 
data are large and so it is certainly possible that the 
measurements by these different methods agree, and 
agree with the theoretical result N^{T BBN ). 

On the other hand these results indicate that 
there is an increase in N V (T) between the end of 
the BBN epoch and recombination. We will show 
that one would naturally expect an increase in N v 
when the effective neutrino temperature approaches 
the neutrino mass. To obtain this result we employ 
that fact that once kinetic (commonly called ther- 
mal) freeze-out has occurred, neutrinos in essence 
free-stream (up to additional e.g. e + e _ annihilation- 
feeding [9]). We must therefore use solutions of 
the Boltzmann-Einstein-Vlasov (commonly called 
Einstein- Vlasov) non-equilibrium kinetic evolution 
equation to describe how the neutrino momentum 
distribution changes in an expanding universe af- 
ter kinetic freeze-out. Such solutions will preserve 
the comoving entropy and neutrino particle number, 
which a Fermi-Dirac integrated equilibrium distribu- 
tion in general does not do. 

Interestingly, an increase in N v constitutes an 
independent measurement of the mass of neutrino (s) 
near to the recombination epoch, which has gone un- 
noticed. The essential physics ingredients for neu- 
trino mass to make a mark are 

1. At least one neutrino has a mass near the 
value of the temperature when the universe be- 
comes transparent to photons at T r = 0.253eV 
e.g. where neutrino density becomes a part of 
WMAP acoustical oscillation analysis [111 ]; 

2. After kinetic freeze-out, but in particular near 
to the T r -epoch, neutrinos free stream through 
an expanding universe. Their momentum dis- 
tribution is thus not in equilibrium. It will be 
seen that the thermal momentum shape can 
only persist when the free-streaming particle 
is massless, as is the case for photons. 

3. The effective neutrino number at any epoch 
where they are already free-stream, e.g. at 
BBN, is known. This allows one to fix the 



magnitude of enhancement factor. Deviation 
from N v — 3 at BBN is not necessary, but 
it affects the predicted value of the neutrino 
mass. 

We stress that the nonzero mass of the neutrino 
combined with the non-equilibrium form of the free- 
streaming distribution offers a novel mechanism for 
increasing the value of N v . 

In section [TT] we discuss the non-equilibrium 
properties of neutrinos. In subsection III Al we de- 
rive the form of the neutrino distribution using the 
Einstein- Vlasov equation. In subsection IIIBI we 
compute various moments of the distribution. In 
subsection III CI we compute the photon reheating 
temperature. In section IIIII we combine our results 
and demonstrate that the effective number of neu- 
trinos rises between BBN and ion-electron recombi- 
nation. In subsection IIII Al we describe the change 
in the effective number of neutrinos while in subsec- 
tion IIII Bl we show that if one incorrectly assumes a 
kinetic (thermal) equilibrium Fermi neutrino distri- 
bution, one finds that N u decreases. We present our 
conclusions and discussion in section IIVI 



II. NONEQUILIBRIUM NEUTRINOS 

A. Einstein- Vlasov Equation in FRW 
Spacetime 

Prior to the neutrino freeze-out temperature, 
Tf ~ 2.3 MeV, collisions are significant and keep 
neutrinos in kinetic (thermal) equilibrium, implying 
that the distribution function of each neutrino flavor 
has the Ferm-Dirac form 

f FD (t,E) = |^ [exp(£/T) + If 1 for T>T f (1) 

where g v = 2 is the degeneracy of each neutrino 
flavor. When the temperature drops below Tf, in- 
teractions no longer occur rapidly enough to keep 
the distribution in thermal equilibrium and the neu- 
trinos begin to freely stream. The neutrinos are no 
longer in equilibrium, and hence kinetic theory must 
be utilized to describe the evolution of their distri- 
bution function. 

The general relativistic Boltzmann equation de- 
scribes the dynamics of a gas of particles that travel 
freely in between point interactions in an arbitrary 
spacetime [l3l - [T6j 

p a d xa f-T^p v d pi f = C[f}. (2) 
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Here V^ v is the affine connection (Christoffel sym- 
bol), / is a function on the mass shell 



(3) 



hence Greek indices are summed from to 3 whereas 
j is only summed from 1 to 3. When collisions 
are negligible C[f] — and all particles move on 
geodesies, yielding the Einstein- Vlasov equation. 

We now specialize to collision free homogeneous 
isotropic cosmological solutions and therefore as- 
sume the flat FRW ansatz for the spacetime metric 



g = dt 2 - a(t) 2 (dx 2 + dy 2 + dz 2 ). 



(4) 



Due to homogeneity and isotropy, the distribution 
function depends on t and p° = E only. Therefore 
Eq. ^ becomes 



Edtf + (m 2 - E 2 )— d E f = 0. 

a 

The general solution to Eq. §5§ is known [l5j |: 



d t a , 



f(t,E)=K(x), x 



D 2 



(E 2 



(5) 



(6) 



where K is an arbitrary smooth function and D is an 
arbitrary constant with units of mass. To continue 
the evolution beyond the freeze out time, tf, we must 
choose K to match at tf the equilibrium distribution 
Eq. CE). 

As mentioned in the introduction, after chemi- 
cal freeze-out there are additional processes that can 
feed into the neutrino distribution. To account for 
this we include a fugacity factor, T, that allows for 
a possible overabundance of neutrinos. The effect of 
T = e a is modeled after that of chemical potential 
\x, except that a is equal for particles and antiparti- 
cles, and not opposite. This means a > (T > 1) 
increases the density of both particles and antipar- 
ticles, rather than increasing one and decreasing the 
other as is common when the chemical potential 
is associated with conservation laws such as lepton 
number. The fact that a is not opposite for parti- 
cles and antiparticlcs reflects the fact that both the 
number of particles and the number of antiparticles 
are conserved after chemical freeze-out, and not just 
their difference. The equality reflects the fact that 
any process that modifies the distribution would af- 
fect both particle and antiparticle distributions in 
the same fashion. 

The use of T to account for the processes that 
feed into neutrinos is exact in the temperature in- 
terval after chemical and before kinetic freeze-out, 
since scattering processes re-equilibrate the momen- 
tum distribution to this shape in order to maximize 



the entropy content. However, it is an approxima- 
tion when the additional particle feeding occurs be- 
low kinetic freeze-out: when using the T-modeled 
distribution, we are making an assumption about 
how such feeding affects the momentum distribution. 
The computed effect of e + e" annihilation- feeding Q 
is small, below 5%, so the potential error of this as- 
sumption is not significant. 
With this in mind, we let 



K{x) 



9u 



1 



8tt 3 T -i e y/x+ m yrf + 1 
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and D = Tfa(tf) to match the Fermi-Dirac 
(FD) distribution at freeze-out. The Fermi-Dirac- 
Einstein- Vlasov (FDEV) distribution function after 
freeze-out is then 



f(t,E) = 



1 



87T 3 T -l e ^(B 2 -m2)/T„HmV^ + -y 



where 



T v {t) 



T f a(tf) 
a(t) 



(8) 



(9) 



Eq. © provides the distribution function that de- 
scribes a gas of neutrinos that have been free stream- 
ing in an expanding universe since they froze out at 
T v (tf) — Tf, with a possible overabundance repre- 
sented by the fugacity T. We will call T v in Eq. (|9]) 
the neutrino background temperature, even though 
the distribution of free streaming particles has a 
thermal shape only for m = 0. This language is, 
however, reasonable since apart from the reheating 
factor of photons due to e + e~ annihilation, which we 
discuss in subsection III C[ T„ is equal to the photon 
background temperature. 



B. Moments of the FDEV Distribution 

Here we compute the stress energy tensor, num- 
ber current, and entropy current associated with the 
distribution Eq. ([8]) 




(10) 

(11) 
(12) 
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fc = /ln(/) + (l-/)m(l-/), / = /. 

9v 
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We first work with the general form of / given 
in Eq. (|6]) and later specialize to the explicit form 
Eq.©. 

Isotropy of the metric and of / in momentum 
space implies that the off diagonal elements of the 
stress energy tensor and spacial components of the 
particle number and entropy currents vanish and 
that the pressure is isotropic. Hence we must com- 
pute 



neutrino flavor 



T w = a A / fEd A p 



(13) 



3 J J E 

n° = a 3 J fd 3 p, 

s o = _ a 3_£^ [ hd 3 p 



I f^d s p, /=l...:{. (14) 

(15) 
(16) 
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where \p\ is the Euclidean norm of the spacial com- 
ponents of p^ and E = p° is given by 

to 2 = E 2 -a(t) 2 \p\ 2 . (17) 

Computing T 00 we find 

POO 

T°° = 4ira 3 / K((E 2 - m 2 )/T 2 )E\p\ 2 d\p\ 
Jo 

I* oo 

= 4na 3 / K(a 2 p 2 /T 2 )(m 2 + a 2 p 2 )^ 2 \p\ 2 d\p\ 
Jo 

[>00 

= 47r/ K(z 2 /T 2 )(m 2 + z 2 )^ 2 z 2 dz (18) 



where we made a change of variables z — a(t)\p\. 
Note that z is the physically measured momentum. 
Similarly 

a roo 

7* = jjj K(z 2 /T 2 )(m 2 + z 2 )-V 2 z*dz, (19) 
n° = 4?r / K(z 2 /T 2 )z 2 dz, (20) 



„o _ 9u 



^ ^ H(z 2 /T 2 )z 2 dz, 



(21) 
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H = K]xiK + (l-K)\n(l-K), K = —K. 

(22) 

We now rename z to p, so that p represents the 
magnitude of the physical momentum, drop the su- 
perscripts, and insert Eq. ((TJ) for K, giving the en- 
ergy density, pressure, and number density for each 



(m 2 + p 2 ) p 2 dp 
2tt 2 7 T -i eV / F r ^(T„ 2 /T|)/T l , + \ : 

(to 2 + p 2 y 1/2 p 4 dp 



P = 



6tt 2 Jo r-i e \ / P 2 + m2 ( T S/ T f)/ T - + i ' 

9u_ r P 2 d P 

27T 2 Jo T -l e yJp 2 +raHT?JTj)/T v + j ' 



(23) 

(24) 
(25) 



These differ from the corresponding expressions for 
an equilibrium distribution in Minkowski space by 
the T v (t) /Tf factor multiplying to only in the expo- 
nential. 

By making a change of variables u — p/T u , one 
sees that both n and s are proportional to T 3 . By 
definition, T„ is inversely proportional to a, hence 

o?n — constant and a 3 s = constant. (26) 

This proves that the particle number and entropy 
in a comoving volume are conserved. This result 
does not depend on the particular form of K chosen 
above. Note further that since an eV scale or be- 
low neutrino mass is at least 6 orders of magnitude 
smaller than Tf, we can ignore the neutrino mass in 
the exponent. This remark does not extend to the 
energy factor multiplying the FD-distribution in the 
calculation of energy density or pressure and there- 
fore Eq. (|2"3")l and Eq. (IM1) have unusual properties, 
while Eq. (|25[) is just what one may naively expect, 
remembering that the mass term in this expression 
is completely negligible. 



C. Photon and neutrino temperatures 

To compare our results with observation, we 
must understand the photon to neutrino tempera- 
ture ratio. The difference between the two is the re- 
sult of the annihilation of e + e~ , mainly into photons. 
Before this reheating, neutrinos, photons, electrons, 
and positrons have the same temperature. When the 
temperature approaches and drops below the elec- 
tron mass, the electrons and positrons annihilate. 
The resulting ratio of photon to neutrino tempera- 
tures has often been studied but not in the type of 
non-equilibrium model presented here. After pho- 
ton reheating, both the neutrino and photon tem- 
peratures evolve inversely proportional to a, so their 
ratio after reheating equals their ratio today. We use 
subscripts 1 and 2 to denote quantities before and 
after reheating, respectively. 
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Before reheating, the entropy in a given volume, 
V\, is the sum of relativistic neutrinos, electrons, 
positrons, and photons 



Si 



o 7 7 

3o9v + o-9e± + .9 7 



2tt2 
45 



(27) 



where 7\ is the common neutrino, e + e~, and 7 tem- 
perature. The particle-antiparticle and spin-helicity 
statistical factors are g v = 2,g e ± — 4, <? 7 = 2. 

As computed in Ref.lgj, not all products of 
e + e~-annihilation go into the photon gas, some go 
into overpopulating the neutrino distribution. We 
model this by an T > 1 and assume that the neu- 
trino temperature is not reheated. After e + e~- 
annihilation, the entropy in the corresponding vol- 
ume is the sum of relativistic neutrinos with T > 1 
and photons 



2^_ 2 
45 



S 2 =['SS, + ^ r g 1 T^)V 2 , 



/"DO 

S v = -^[T»,2j H{u 2 )u 2 du 



(28) 



(29) 



where H is given by Eq. (J22J) . Note that we must 
now distinguish between the neutrino and photon 
temperatures. The assumption that the neutrino 
temperature is not reheated implies that 



Ti/T„ i2 = oa/ai = (1-2 M) 1/3 ■ 



(30) 



Equating the two entropies and using this fact we 
obtain 



1/3 



(31) 



R = 21 135 r H{u 2 )u i du _ 

8 47T 4 i 



Writing T = e CT and expanding this to first order 
about a = gives 



1/3 



T. 



7,2 



1 



405C(3) 
22tt 4 ' 



(32) 



where is the Riemann zeta function. However, 
Eq. (|32|l is only applicable for relatively small a. Fig- 
ure [T] shows the ratio Eq. (|3"Tj) as a function of T as 
obtained numerically, where we have ignored terms 
involving m v /T f ~ 10~ 6 < 1. Note that for T = 1, 
we obtain the usual value T lt2 /T Ut2 = (11/4) 1/3 . 
Equal photon and neutrino temperature is obtained 
from T = 2.11. The maximum value of T, for 
which all entropy goes into neutrinos, is approxi- 
mately T = 2.91. 




FIG. 1: Ratio of photon to neutrino temperature after 
reheating as a function of T. 



one massive neutrino 
two massive neutrinos 
three massive neutrinos - 




FIG. 2: Maximum neutrino mass as a function of T, 
obtained by equating neutrino and dark matter densities. 



Using the measured CMB temperature and the 
above photon reheating calculation we obtain the 
current neutrino background temperature as a func- 
tion of T. In turn, this allows us to set a bound 
on the mass of the heaviest neutrino(s) as a func- 
tion of T by equating the neutrino energy den- 
sity, Eq. (|2"3"1) . to the measured dark matter density, 
Pdm = 8.934meV 4 [T3. We show the results in 
the case of one, two, or three equal mass neutrinos, 
where we set the mass of the lighter neutrino (s) to 
zero in Figure O The maximum allowed mass goes 
to zero in Figure [2] near to the maximum value of 
T from Figure [1] for a large range of values of pdm- 
This is because for fixed T 7 , T v diverges to infinity 
as T approaches its maximum value, and hence the 



G 



neutrino energy density diverges to infinity as well. 
As we will see, we do not violate the dark matter 
limit in any of these three cases, since the value of 
T that is of most interest to us is near to T ~ 1.1, 
and the neutrino masses we consider are near to or 
below 1 eV. 



III. EFFECTIVE NUMBER OF 
NEUTRINOS 

A. Non-equilibrium model 

Typically, when considering the evolution of the 
Universe the neutrino distribution is modeled as 
massless. The difference between the massless model 
and the more precise model based on the Vlasov 
equation that we derived above can be quantified 
as an effective number of neutrinos, defined by com- 
paring the total energy density or pressure to a gas 
of massless particles and antiparticles 

E^= ^ 2N »AT„V (^) 4/3 T 7 4 , (33) 

i 

5> 55 3lT 2A ^ (r - T) (n) 473 ^' (34) 

where the sum is over the three neutrino flavors and 
(4/11) 1 / 3 is the reheating factor when T = 1, as is 
normally assumed. There are two different versions, 
depending on whether we define N v via pressure or 
energy density. The value given in Q is defined using 
the neutrino energy density, and so we will do the 
same. From here on, N v will refer to N VtP . 

To compare the energy density Eq. ([23")) to that 
of a massless particle distribution, we make a change 
of variables u = p/T v and neglect terms involving 
m/T f < 1 



2tt 2 



(m 2 /T 2 + v? 



,1/2 2 



T 1 exp(u) + 1 



-du, (35) 



where the upper index 'EV reminds us that we have 
used Einstein- Vlasov free streaming solution for the 
neutrino distribution. For m <C T„ this approxi- 
mates the energy density of the equilibrium distri- 
bution of a massless fermion, but when T v is on the 
order of the mass, the mass term becomes important 
and changes the effective number of neutrinos. This 
is in addition to the increase in effective number of 
neutrinos due to T ^ 1. 

To explore the small mjT v limit, let T = e a 
and a = m/T„. As an illustration of the qualitative 



behavior, we Taylor expand Eq. (|3"5j) . Note that the 
following expansions will not be valid to the required 
precision in the regime where we will need numerical 
results 



-ir 2 T? 
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n 2 cr 2 
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where C(3) = J2T=i 1 / k 
approximation regime 
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4/3 
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1.20206. 
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(36) 
In this same 
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where N — 3 is the number of neutrino flavors in 
the standard model and the sum is over the three 
neutrino species. From this expansion, we can see 
the claimed increase in N„ as the temperature de- 
creases. We now give a more explicit approximation 
that is useful in the region of parameter space that 
we are most interested in. The coefficients of the 
quantities involving only a are taken from Eq. (|37p , 
as that expansion is valid for the values of a of inter- 
est when a = 0. This allows us to match the large 
T 7 asymptotic behavior with a high degree of accu- 
racy. The remainder of the coefficients are obtained 
via a nonlinear least squares fit to N^ v /N 

N EV fh 2 
=1 + 1.861cr + 2.298<7 2 + .0296— r 
N T 2 



2.493cr 3 + .0357cr^- + 2.574cr 4 

1 

7 

.0618^^-4.7910-^^, (38) 



T 2 

7 



rp4 

7 



where m — m i * s the neutrino mass radius. 

The fit is valid to within 3.6% relative error in the 
region 1 < T < 1.2, < mi/T 7 < 5. The coefficients 
differ from those in Eq. (|37l) since we have focused 
on achieving a fit near m/T 7 ~ 1 where the effect we 
are investigating is most pronounced. In addition in 
Eq. psp , we have converted from neutrino to photon 
temperature using Eq. pip to facilitate comparison 
with measured values. 

The effect of mass becomes important when the 
neutrino temperature approaches the mass of the 
heaviest neutrino. For illustrative purposes, a plot 
of N u /N for three massive neutrino of mass .35eV 
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T T [eV] m [eV] 



FIG. 3: Effective number of neutrinos derived from the 
energy density for three neutrinos of mass .35eV with 
fugacities T from 1 to 1.2 (solid lines bottom to top in 
increments of .05). The series expansions (doted lines) 
and nonlinear fit (dashed lines) are also shown. The 
vertical line marks the photon temperature at recombi- 
nation T r = .253eV. 



FIG. 4: Comparison of N v at recombination to N v at 
BBN for (bottom to top) one, two, and three massive 
neutrinos with fugacity T = 1.12. The horizontal lines 
show the values corresponding to N BBN = 3.71 (bottom 
line) and Nt hBBN = 3.046 (top line). 



fit for the coefficients 



N, 



r.EV 



(solid lines) and various values of T is given in Figure 
[31 We also show the series expansion Eq. (f3"7f (dots) 
and the nonlinear fit Eq. (|38p (dashed lines) . These 
graphs are approximately constant for T 7 3> T r , and 
so the values can be extrapolated back to BBN and 
are independent of the neutrino masses. From the 
high temperature expansion in Eq. (|37l) we see that 
to match 7V„ = 3.71 at BBN we require a fugacity 
of T = e a = 1.12. If instead, the value computed in 
Q is correct and the discrepancy is simply measure- 
ment error, we require a fugacity of T = 1.01. 

We can also expand the ratio of effective num- 
ber of neutrinos at recombination to that at BBN, 
N^/N BBN , as a series in a and a. To third order 
this equals 



N, 



r,EV 



jyBBN 



1+ - 



2/3 — 2 
J^2 



21tt 2 
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where we have used Eq. (|3lj) to replace the neutrino 
temperature with the photon temperature at recom- 
bination. We again perform a nonlinear least squares 



N BBN 



.0281 



rp2 



.00122 er 



JT 2 



_ 3 .6810- 4 a 2 ^-3.9910- 4 ]T^. (40) 

7 i 7 

The fit is valid to within 3.0% relative error in the 
region 1 < T < 1.2, < TOj/T 7 < 5. Figure 0] shows 
the variation of the ratio N^/N BBN as a function 
of the neutrino mass, in the cases of one, two, or 
three equal mass neutrino(s), setting the mass of the 
lighter neutrino(s) to zero. The result obtained by 
computing the ratio numerically using Eq. pip and 
Eq. (f33f is shown in solid lines and the nonlinear fit 
is shown in dashed lines. Note that at higher mass, 
we see the power law in to start to break down as 
higher order terms in the expansion become impor- 
tant. The result is shown for T = 1.12, but it is 
relatively insensitive to T. This can be seen from 
Eq. (f3"9"]> . since at lowest order there is no dependence 
on a 

For the three possible cases of Nm = 1,2,3 (col- 
umn 1 in table H|, the values of the neutrino mass 
m v , which result in a rise of from ^ hBBN — 
3.046 (i.e. T = 1.01) or N BBN = 3.71±°;g (i.e. 
T = 1.12±g;{$) to N r u = 4.34l ^ are given in col- 
umn 4. Depending on the data and analysis used, 
the neutrino mass values we find are close to, or out- 
side of current cosmological bounds on the neutrino 
mass from WMAP data 181 ■ However, these bounds 
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TABLE I: Nm is the number of massive neutrinos as- 
sumed. The neutrino mass, m„ is determined by con- 
necting N v (Tbbn) to N v {T r ). The last column shows 
the present day (PD) contribution of neutrinos to dark 
matter energy density, with T and the neutrino masses 
fixed to describe the increase of neutrino degrees of free- 
dom from BBN to recombination epoch. Notation b c a 
means that c is the central value bracketed by (o, b). 



N M 


N v {Tbbn) 


N v (T r ) 


m v [eV] 


PD / PD 
Pv 1 Pdark 


1 


3.046 


A 34+-086 

^-^-o.ss 


1,66 1.13o.528 


° 16 0.11o.05 


2 


3.046 


A Q 4 + .086 
^■•3^-0.88 


°-" 5 0.703o.347 


19 0.14 .o7 


3 


3.046 


4 Q4+-086 
^■•3^-0.88 


°' 754 0.542 .275 


°' 22 0.16 .o8 


1 


o 71 +0.47 
' i -0.45 




1 ' 47 0.627o.o 


15 0.07 210 -4 


2 


n 71 +0.47 
' 1 -0.45 


4 Q4+-086 


°' 892 0.408 .o 


18 0.09 210 -4 


3 


o 71 +0.47 
°- ' ± -0.45 


4 34+ 086 

^■■5^-0.88 


°- 680 0.322 .o 


°' 21 0.10 210 -4 



could change based on the effect of neutrino mass 
and non-equilibrium distribution described here. A 
more robust bound for our purposes is the non- 
cosmological bound m v < 2.05eV [ljj which we do 
satisfy. 

We note that the scale of m v is non-negligible 
compared to the limit set in Figure[2j where we asked 
how large a mass is needed to describe all dark mat- 
ter in the present day (PD) Universe. As a conse- 
quence, if the increase in N v central values were to be 
confirmed then the neutrinos emerge as a significant 
component in the dark matter inventory, alte ring the 
presumed energy inventory in the Universe [19j. In 
the last column in table fl] we present the fractional 
contribution of neutrinos to the present day dark 
matter component for all cases. Note that dark mat- 
ter is 22% of the present day (PD) energy density in 
the universe (l7j . 



B. Problems with Equilibrium Distribution 

Here we show that a naive application of equi- 
librium thermodynamics will produce a reduction of 
N v and remark on other issues that arise in such 
an approach. If one assumes an equilibrium dis- 
tribution of massless neutrinos, or more precisely, 
negligible neutrino mass as compared to the ambi- 
ent temperature, then the neutrino energy density 
scales as Ti, and hence the effective number of neu- 



trinos, as defined in Eq. (|33[) . is constant in time. 
Now suppose we use the fact that neutrinos are mas- 
sive, but neglect the non-equilibrium FDEV aspects 
in the momentum distribution of particles. To make 



this mistake, consider the energy density of each neu- 
trino flavor to be, instead of Eq. (|35[). 



P 



Pi — 



9^2 
2vr 2 



i 2 ) 1/2 u 2 du 



o T _1 exp(V a 2 + u 2 ) + 1 



(41) 



Expanding this about a = = m and a = In T = 
to third order gives 



J_ 

120' 



7T J„ 



120 (<b „. . ir 2 a 2 
1+ 7^ 2 C(3)a 



24 



ln2<7 3 - — ^j— Q^(7 



(42) 



and hence 
N c<1 / 11 \ 4 / 3 j 14 



N 



rp4 

7 



E 



40 /9 



7tt 4 V 2 
In 2 



C(3)a 



24 



+ ln2er 3 a 2 a 

2 1 



(43) 



Comparing this with the expression derived from the 
non-equilibrium distribution, Eq. (|37|) . we find the 
opposite qualitative behavior in a. 

Thus basing the analysis on a naive equilibrium 
distribution, one would predict that N v decreases 
with temperature rather than increases. Hence 
one arrives at the incorrect qualitative behavior of 
N v unless both the neutrino mass and the non- 
equilibrium form of the neutrino distribution are ac- 
counted for. 

The reader should note that there are other 
immediate problems with an equilibrium approach. 
The FD distribution does not conserve with dimin- 
ishing temperature the comoving entropy and parti- 
cle number content when ra/0. Even if we manipu- 
lated the value of T to preserve the particle number, 
we cannot also preserve the neutrino entropy since 
it is the nature of thermal equilibrium that entropy 
is exchanged with the heat bath and we would need 
to manipulate T to be able to preserve entropy, but 
that is controlled by other forces. We conclude that 
use of FD equilibrium distribution has multiple in- 
consistencies and in its most naive application leads 
one to predict a reduction of N v in an expanding 
Universe. 



IV. DISCUSSION 

After the freeze-out, neutrinos freely stream 
through the expanding universe and their distribu- 
tion function ceases to be thermal. Using a solution 
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of the Einstein- Vlasov equation, we find that as the 
universe cools and the temperature approaches the 
mass of the heaviest neutrino, the effective number 
of neutrino degrees of freedom, N u , increases. Such 
an increase is favored by the measurements of N v 
at BBN and electron-ion recombination. We have 
shown that this increase is due to the interplay of 
the neutrino mass and the non-thermal form of the 
neutrino distribution after freeze-out. Both these 
properties are essential for reaching the correct con- 
clusion, i.e. that there is an increase in the effective 
number of neutrino degrees of freedom. 

It is helpful to understand the essential physics 
of this observation without the need for much the- 
ory. In order to move from the number density, 
Eq. (|25l) , to the energy density, Eq. (l23t , we weight 
the energy of a neutrino E u = \/p 2 + m 2 by the 
non-equilibrium free streaming FDEV distribution 
Eq. (JSJ . The distribution preserves both the comov- 
ing particle number and entropy and, in case of neu- 
trinos, the mass term in the distribution can always 
be neglected. However, one cannot neglect the mass 
in the energy of a neutrino E v , as at some point the 
momentum of a neutrino in an expanding universe 
will be small. Thus one weights a massive parti- 
cle energy with an effectively massless distribution 
function at a temperature T v . It is possible to ar- 
rive at this result with entirely intuitive reasoning, 
the exact theoretical derivation is the essence of this 
report. 

The values of m„ seen in table fl] are close to 
or exceed various cosmological bounds [lq . How- 
ever, these bounds could change when our non- 
equilibrium FDEV form Eq. (|8) of the free-streaming 
neutrino distribution is incorporated in the analysis. 
Recall that there is no direct laboratory measure- 
ment of the magnitude of the neutrino mass except 
for the limit on electron-neutrino m„ < 2 eV derived 
from Tritium decay [TtJ which we satisfy. The differ- 
ences between neutrino masses (squared) have been 
derived from observations of neutrino mixing, i.e. 
Am 2 12 = 75 ±2 meV 2 , Am| 3 = 2320^0° meV 2 
Therefore some increase in N v between BBN and 
recombination is unavoidable. 

Due to the large uncertainties in the measured 
N v values, they are still (barely) consistent with a 
scenario with three neutrinos whose mass is neg- 
ligible, relative to the recombination temperature 
(ra v <C T r ), and therefore results in a negligible in- 
crease in N v . Even so, note that to be consistent 
with the given error bounds, this requires 3.46 < 
]\[Bbn < 4 2.8 which is larger than the calculated 
value and corresponds to 1.08 < T < 1.19. Thus, an 
important observation is that if one wants to main- 



tain the calculated value N* hBBN ~ 3.046 @, then 
a neutrino mass on the order of the recombination 
temperature is required to account for the increased 
value N£. 

To satisfy the often cited cosmological limit 
J2 m v < 0.6 eV (this could be subject to change 
considering our results), and considering the light- 
est central value in table Q] (m v = 0.331 eV), we see 
that all three neutrinos would need to be around 
m v ~ 0.2 eV in mass (i.e. m 2 = 40,000meV) in or- 
der to easily satisfy the measured values Am 2 . For 
the case of three neutrinos of mass .2eV, we predict 
N^/N BBN ~ 1.07. This is consistent with the mea- 
sured range 0.83 < N^/N BBN < 1.60. In general, 
we conclude that an increase in N v signals that the 
mass of the heaviest neutrino is near the recombina- 
tion temperature T r . 

To conclude: We have shown that the observed 
increase in effective number of neutrinos as the tem- 
perature decreases is a natural consequence of the 
presence of non-equilibrium, free streaming FDEV 
neutrino distribution Eq. (|5J) after freeze-out. The 
deviation of the energy density of the FDEV non- 
equilibrium distribution from that of a Fermi-Dirac 
equilibrium distribution is due to the finite mass 
of the neutrinos. Within the standard model, the 
increase in effective number of neutrinos N v over 
time leading to the recombination epoch depends 
solely on the neutrino mass and presence of non- 
equilibrium distribution. 

In table U we have presented explicit values for 
m v favored by the current knowledge of the values 
of N v . Note that the smaller the number of domi- 
nant mass neutrinos and the larger the change in N v 
needed between BBN and recombination, the larger 
is the value of m v we find. If little change in N v 
is needed, m v /T r — > 0. However, within the cur- 
rent error bounds, it is not possible to have negli- 
gible neutrino masses (m„ <C T r ) when using the 
calculated value of N t J lBBN Q, see the first three 
lines in table HI where each entry has a lower finite 
m„. If the central values of N u {T r ) prevail within 
error bars, and no new mechanism is discovered to 
further increase N v (Tbbn) then the relic neutrinos 
will have a mass above 0.275 eV and will contribute 
between 5% and 22% of the dark matter inventory 
in the Universe, or if relic neutrinos are counted as 
visible matter, they will increase the visible matter 
inventory noticeably, with the potential to double it. 
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